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ABSTRACT 

We present a multi-frequency far-field beam map for the 5m dish telescope at the Bleien Obser¬ 
vatory measured using a commercially available drone. We describe the hexacopter drone used in 
this experiment, the design of the flight pattern, and the data analysis scheme. This is the hrst 
application of this calibration method to a single dish radio telescope in the far-held. The high 
signal-to-noise data allows us to characterise the beam pattern with high accuracy out to at least the 
4th side-lobe. The resulting 2D beam pattern is compared with that derived from a more traditional 
calibration approach using an astronomical calibration source. We discuss the advantages of this 
method compared to other beam calibration methods. Our results show that this drone-based 
technique is very promising for ongoing and future radio experiments, where the knowledge of the 
beam pattern is key to obtaining high-accuracy cosmological and astronomical measurements. 

Subject headings: radio, calibration 


1. INTRODUCTION 

In the next decade, a number of large radio experi¬ 
ments are scheduled to begin data collection. One of 
the key science goals of these programmes is to map 
the HI intensity in the Universe through its 21 cm emis¬ 
sion line. HI intensity mapping provides a probe of the 
Baryon Acoustic Oscillation (BAO) feature in the matter 
power spectrum (or more directly, the HI power spec¬ 
trum) that is independent of traditional measurements 


intensity mapping are still dominated by ~ 20% statisti- 


cal errors ( 

Padmanabhan et al. 

2015 

), with future large 

surveys where percent-level sta 

istical uncertainties are 


expected, the systematic errors from the knowledge of 
the telescope beam needs to be controlled to sub-percent 
level. 

Traditionally, beam calibration has been done using 


brigh t astronomi cal sources such as the sun (e.g. Kraus 


19661, the moon (Tello et al. 20131, and known bright ra- 


using galaxy clustering and weak leasing (Wyithe et al. y-, 

2^ IChang et al.pOlOl iBulletal. 1120 151) Examples of ^nd Virgo A ( |Baars et al |1977[ ). Having the source drift- 
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dio sources s uch as Cassiopeia A, Taurus A, Cygnus A, 


Sky Survey (HIPASS, 

Zwaan et al. 2005), 

the HI Jo- 

drell All-Sky Survey (1 

Hl-JASS 

Lang et al. 

\ / A J J 

2003 

\ 1 1 W 

), and 


tin et al.]|2010). Many more future programmes are also 


being denned and built, including the Square Kilome¬ 
tre Arra }[T (SKA), the Low Frequency Arraj0 (LOFAR), 
and the Baryon aco ustic oscillations In Neutr al Gas Ob¬ 
servations (BINGO, Battye et al.||2012] |2013| ). 

For HI intensity mapping, especially at low redshift, 
an advantageous survey configuration is to operate dish 
arrays of small to moderate dish sizes (5-15m) in single¬ 
dish configurations (jSantos et al.|2015|. This allows wide 
collecting area and a complete sampling of relevant spa¬ 
tial scales. However, in order to achieve the required 
accuracy for the single-dish telescopes, one needs to un¬ 
derstand and calibrate the response pattern, or the beam 
of each telescope very well. Small deviations of the me¬ 
chanical configuration or the environment (e.g. temper¬ 
ature, wind) can cause changes to the beam pattern and 
introduce systematic errors in the measurement. A 1% 
uncertainty in the size of the beam roughly propagates 
into a 4% systematic uncertainty in the amplitude of the 
power spectrum. Although current constraints from HI 

* Electronic address: chihway.chang@phys.ethz.ch 
^ https://www.skatelescope.org/ 

^ http: / /WWW . lof ar. org 


scan over the extent ot the beam gives one a measure of 
the beam shape convolved with the source. Similarly, 
one can use satellites or other artificial sources placed on 
distant towers to perform such calibration. However, in 
the case of astronomical sources, the number of usable 
sources is limited and decreases for smaller radio dishes. 
Furthermore, the fluxes and the sizes of these sources 
can fluctuate over time. In the case of satellites, the 
frequency range of the source spectrum is usually very 
limited, though the intensity is fairly high and regular. 
To avoid these limitations in using astronomical objects 
or satellites as calibration sources, the ideal solution is 
to construct a artihcial calibration source that is flexi¬ 
ble and controllable so one can tailor it to the specific 
telescope and experiment of interest. 

In this paper, we implement this idea by using a noise 
source carried by a commercial hexacoptor drone. The 
noise source emits a flat spectrum in the frequency range 
980 MHz - 1250 MHz at high power. The frequency 
range is chosen to be the 21 cm frequency redshifted to 
z = 0.14-0.46. The drone flies in a region where the far- 
held beam pattern can be mapped. We show that this 
method gives a controllable, light-weight solution to the 
beam calibration problem for radio telescopes. An earlier 
stud}0with a similar setup has been done by the aper- 

® http://www.skatelescope.org/wp-content/uploads/2011/ 
03/SKA_NEWSLETTER_V0LUME_25.pdf 
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Fig. 1.— Schematics of the experiment setup. The drone, indi¬ 
cated by the green circle, is flying about 150 m above ground, in 
a plane directly above the telescope. Note that the dimensions in 
the plot are not drawn to-scale. 

ture array verification programme (AAVP) as a proof of 
principle. In this paper we work with a different telescope 
type and wavelength. We also perform more quantitative 
analyses of the data to show the potential of this method. 
Note that although this work is motivated by HI inten¬ 
sity mapping cosmology, the application of our method 
can be extended to other science areas where single-dish 
radio telescopes are used in the centimeter wavelength, 
such as solar physics, pulsars and radio bursts. 

This paper is organised as the follows. In ^ we de¬ 
scribe the separate components of the experiment de¬ 
sign. In ^we introduce definitions of the characteristic 
quantities we like to measure from our beam. The data 
processing, analysis and hnal results of this experiment 
are presented in 0 We also compare our measurements 
from the drone with other more traditional approaches. 
In ^ we present our conclusions. 

2. EXPERIMENT DESIGN 

We describe below the separate components of the ex¬ 
periment: the telescope, the spectrometer, the drone, 
the noise transmitter, and the design of the flight pat¬ 
tern. Figure [J illustrates the schematics of the experi¬ 
ment setup. A test run of smaller scale was performed 
on October 28, 2014, while the full experiment was car¬ 
ried out on November 20, 2014. All data presented in 
this paper are from the latter dataset. 


2.1. The Bleien 5m dish 


We demonstrate our new calibration technique on the 
5m parabolic dishj^ (f/D= 0.507) at the Bleien Obser- 
vatorjldin Granichen, Switzerland (longitude 8.112215°, 
latitude 47.3412278, altitude 469 m). The surrounding 

^ The radio dish was constructed in 1972 by the Swiss company 
Schweizerische Wagons und Aufziigefabrik AG. 

® http://www.astro.ethz.ch/research/Facilities/ 

Radioteleskop_Bleien 


Fig. 2.— Image of the 5m radio dish at the Bleien Observatory. 
From the image one can see the wire meshed reflector, as well as 
the four struts which hold the horn receiver in focus. 


area ^ 1.5 km in radius is protected against commercial 
radio emission in the frequency range of interest for this 
experiment. Figure shows an image of the telescope. 
During the experiment, the telescope is pointed vertically 
up towards the zenith so that the beam pattern can be 
me asured at a plane parallel to the ground. As discussed 
in §4.6.2[ the telescope can also be pointed at different 
elevation angles to map the beam shape as a function 
of elevation. A cylindrical horn feed is supported via 
four struts at the focal point. The cylindrical horn has 
a length of 185 mm, diameter of 200 mm and the dipole 
length is 85 mm. 


2.2. The CALLISTO spectrometer 

During the experime nt, data is collected by the CAL¬ 
LISTO spectrometer (jBenz et al. 2005). The CAL¬ 
LISTO spectrometer is a programniable heterodyne re¬ 
ceiver built in the framework of IHY2007 and ISWI by 
former Radio and Plasma Physics Group at ETH Zurich, 
Switzerland. The instrument natively operates between 
45 and 870 MHz and has a frequency step size of 62.5 
kHz. The data from the telescope is down-converted to 
match the frequency range of CALLISTO. The data ob¬ 
tained from CALLISTO are FITS-files with up to 400 
frequencies per sweep. We set the time resolution of the 
data to be 0.25 sec and 200 channels per spectrum. The 
integration time is 1 ms, the radiometric bandwidth is 
^300 kHz (~ 5 times the frequency step size), and the 
overall dynamic range is larger than 50 dB. 


2.3. Drone 


The hexacopter drone used in this experiment was pro¬ 
vided by the private company Koptershoi0 The main 
characteristics of the vehicle are listed in Table while 
an image of it is shown in Figure The critical features 
of the drone considered in this work are the following: 


http://www.koptershop.ch/ 
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Fig. 3.— Image of the drone and noise transmitter horn used in this experiment. The left image is taken during the flight and the right 
image shows a zoom-in of the vehicle, transmitter and gimbal. The red ball on the drone is the “nose”, which helps the user to identify 
the orientation of the vehicle. The noise transmitter horn can be seen held by the gimbal and pointing vertically down at all times. 


TABLE 1 

Basic characteristics of the drone. 


Quantity 

Specification 

Diameter of full vehicle 

110 cm 

Weight 

10.88 kg“ (total) 

Maximum motor power 

2.01 kW 

Propeller dimensions 

16” (diameter) X 6” (pitch) 

Flight control system 

DJI WooKong-M'> 

Maximum flight duration 

13.5 minutes 


Including 2.73 kg for the weight of the accumulator. 
^ http://www.dji.com/product/wookong-m 


• The drone should be able to carry the weight of the 
noise transmitter in addition to its own weight. 

• The gimbal on the drone should be able to steadily 
point the noise transmitter to a given direction, 
which means that the drone flight needs to be sta¬ 
ble, and the gimbal should compensate for any in¬ 
stability. 

• The drone should be able to sustain a flight long 
enough for at least one pass through the expected 
beam pattern. 

• The 3D position of the drone should be recorded to 
sufficient accuracy when the signal is transmitted 
from the noise transmitter. 

All the above requirements can be met with commercially 
available drone vehicles. Specifically, the gimbal used 
was purchased commercially and modified to fit the spec¬ 
ifications of the experiment. For the last point above, the 
position of the drone during the flight is given by 5 or 
more GPS satellites (depending on the situation during 
the flight). The GPS absolute positioning is accurate to 
a few meters in the transverse direction (we estimate the 
relative accuracy in ®. The height of the drone posi¬ 
tion is controlled by a oarometric altimeter carried by the 
drone and is re-calibrated before each flight. The baro¬ 
metric altimeter measurements are typically accurate to 
well within a meter. 

2.4. Noise transmitter 


The noise transmitter is composed of a non-coherent 
semiconductor noise source, an attenuator, a broad-band 
amplifier, a band-pass filter, a power supply and a trans¬ 
mission antenna. The whole unit is light-weight (< 2 
kg), making it possible to be carried by the drone for 
extended flights. 

The band-pass filter ensures that only the frequency 
range of interest is transmitted. Approximately 3W of 
total power is needed for the noise transmitter, which is 
separate from the power supply for the drone. Finally, 
the transmission antenna is a double ri^ed horn antenna 
as shown in the right panel of Figure!^ The antenna is 
constructed with thin light-weight aluminium sheets and 
covered by a polystyrene plane, with gain of maximum 5 
dB and a frequency range similar to our band of interest. 
The antenna is linearly polarised, with a fluxgate mag¬ 
netometer to maintain the stability of the polarisation. 

Note that to transmit in our frequency band at such 
intensity in Switzerland, a transmission permit from the 
Federal Office of Communications (OFCOM) was ob¬ 
tained. 

2.5. Design of the flight pattern 

The flight pattern is designed to fully cover the extent 
of the beam at far-held, while having sufficient resolution 
of the high-order side lobes. We hy at an elevation of 
about 150 m. At the frequency of interest (^ 1 GHz), 
this elevation is sufficiently close to the far-held region 
dehned by the Fraunhofer distance: 

. 202 

df = —— « 166 m, (1) 

A 

where D is the telescope aperture and A is the wave¬ 
length of interest. Equation suggests that for larger 
telescopes and shorter wavelengths, the far-held require¬ 
ment is more stringent. The commercial drone we used 
provide an appropriate platform for testing the Bleien 5m 
telescope at these wavelengths. For experiments with a 
much larger far-held requirement, more advanced drones 
can be used. Alternatively, one can carry out the experi¬ 
ment in near held and reconstruct the far-held beam via 
modelling (an example modelling software is introduced 
in Appendix [A|) . 
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Given the elevation mentioned above, we set the flight 
pattern to be on a rectangular grid of 75 mx 75 m directly 
above the telescope, which comfortably covers the beam 
out to the fourth side lobe. The grid is oriented in the 
North-South (NS)/East-West (EW) direction with each 
flight track separated by 5 m. The 5 m spacing corre¬ 
sponds to about 1.9° at the beam centre 0.5 times the 
FWHM of the beam at IGHz), suggesting roughly three 
tracks would pass through the extent of the main beam. 
This means there are 16 flight tracks in the NS direction 
and 16 in the EW direction. The flight time limit of the 
drone allows it to complete two tracks for each flight, 
then the batteries need to be changed. The polarisa¬ 
tion angle of the reciever feed horn is 45° in the NW-SE 
direction, while the polarisation angle of the noise trans¬ 
mitter is always parallel to the direction of flight. This 
suggests that there should be no amplitude differences 
for the tracks along the NS and EW tracks due to the 
polarisation. The left panel of Figure H illustrates the 
schematics of the raster scan pattern for this experiment. 
The full flight pattern is programmed into the drone con¬ 
trol software so that it runs automatically. Manual con¬ 
trol is invoked only during takeoff and landing for safety 
considerations. 


3. BEAM CHARACTERISATION 

The radio beam prescribes the sensitivity of a radio 
telescope as a function of the angle of the incoming ray 
relative to the telescope pointing. Typical beam profiles 
are composed of a prominent main beam and side lobes. 
The effect of the side-lobes is to pick up signals that 
are not in the direction of interest. The larger the main 
beam is relative to the side lobes the more efficient the 
beam is in collecting signal. The narrower the spatial 
extent of the main beam, the higher its resolution. The 
goal of this paper is to map the 3D (2D in angular space 
-|- frequency) beam of the Bleien 5m radio dish. From 
the beam map, we also derive basic characteristics of the 
telescope. 


We follow the terminology used in Rohlfs (1986). As¬ 


sume P{0', A) to be the normalised beam with peak in¬ 
tensity equal to one and falls to zero at infinity. The 
first convenient measurement is the Full-Width-Half- 
Maximum (FWHM) of the main beam, or the average 
diameter of the contour where P{0; A) = 0.5. In an ideal 
case with only a perfect circular top-hat aperture, the 
FWHM depends on the wavelength A and the dish diam¬ 
eter D according to 


FWHM(A) = 1.028^. 


( 2 ) 


Second, we can integrate P over the full 47r solid angle 
to get the beam solid angle or integrate only inside 
the first null to get the main beam solid angle Hm. 


nA{x) = 


p{0-x)d'^e- 


477 


( 3 ) 


^m{X) = I / P{0-X)d^0. (4) 

/ main lobe 


From flA and we can calculate two other quanti¬ 


ties. The beam efficiency is defined as 


VmIX) 


Hm(A) 


( 5 ) 


dM{X) is a measure of the relative level between the main 
beam and the side lobes. The closer ? 7 m(A) is to 1, the 
more prominent the main beam is and the more efficient 
the beam is in collecting the signal. The effective aper¬ 
ture of a beam is defined as 

Ae(A) = (6) 

The aperture efficiency is defined as the ratio of Ag to 
the geometric aperture Ag = 'k{DI2)'^, or 


??a(A) 


Ae 


4A^ 

tt^aD"^ 


( 7 ) 


In (4.5 we calculate FWHM(A), rjMiX) and ? 7 a(A) for 
our beam measurement. 


4. ANALYSIS AND RESULTS 


In this section we describe the analysis procedure and 
show the results of our beam meas urem ent. We first de¬ 
scribe the data processing steps in §4.1[ Nex t we present 
t he re sults in terms of the ID profile ((4.2), 2D profile 
( ^4.3[ ), and wavelength-dependence ( p.5[ ) of the beam. 
Finally we co mpare the drone measurements with other 
approaches in §4.6[ 

4.1. Data processing 


From the spectrometer, we read out time-series signal 
from the receiver over 200 frequency channels. The first 
task is to match the signal from the spectrometer at ev¬ 
ery time-instant to the drone location in the air when 
this signal was emitted. From the GPS data, we have a 
coordinate record for each emitted signal, which includes 
2-sec “on” and 1-sec “off” signal from the noise transmit¬ 
ter. The right panel of Figure |5 shows the actual GPS 
records for the position of the drone at each pulse from 
the noise transmitter. We calculate the median RMS 
scatter in the longitude (latitude) direction for all NS 
(EW) tracks to be 0.55 (0.51) meters, which corresponds 
to 0.21° (0.2°) at the beam centre, or ~ 5% of the beam. 
All the emitted signals are recorded by the spectrometer, 
suggesting none of the emission signals were too weak to 
be detected. As the drone does not fly at a constant 
speed, the distance between each on-off signal changes. 
We take this into account when assigning a coordinate 
in the air to each of the spectrometer data pixels. 

The raw data from the spectrometer appears as a series 
of on-off signals, with often artefacts at the edge of the 
off transition due to the electronics in the transmitter 
as shown by the green arrows in the top left panel of 
Figure]^ The on-off signal is designed to help remove the 
background and low-level RFI from the signal of interest 
as explained below. The data is cleaned via the following 
steps: 

1. Gonvert the units of the raw data into dl^by mul- 


^ Note that, in this experiment, we can only measure the relative 
level of the beam intensity at each position, as the instrument 
noise and throughput is not calibrated. In the future with better 
characterisation of the noise transmitter, we can consider doing an 
absolute (radiometric) measurement. 
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tiplying the data with the conversion factor used in 
the CALLISTO spectrometer 255(ADU)25?mV/dB) ■ 


2. Match the GPS positions (in longitude and lati¬ 
tude) to the signal received by the spectrometer, 
then convert the longitude and latitude to angles 
from the optical axis of the beam. 


3. For each frequency, subtract all data by the me¬ 
dian value over time. This step removes the time- 
independent l ow-l evel standing-wave pattern (see 
discussion in (4.51. 


4. Remove the “off” signal and the artefacts around 
the “off” signal by first placing a cut at 0.1 dB and 
then removing 4 pixels on each side around the cut. 
Manually mask any remaining “off” signals that 
were not cut oul[3 


5. Linearly interpolate over the gaps in the signal. 

6. Using a similar approach as above, we can remove 
the “on” signal to get the background, including 
the interference from the drone motors. 


7. Subtract the interference from the drone motor 
from the total signal, and rescale the amplitude 
so that the peak is at 0 dB. 

Figure [^illustrates these different steps. 

Two geometrical issues also need to be considered to 
calibrate the reconstructed beam. First, we check that 
the height of the drone during the flight has been stable 
within 10 cm (upper bound of barometric altimeter pre¬ 
cision), this corresponds to a <0.3% change in intensity. 
Second, as the beam pattern is measured on a tangent 
plane at the centre of the beam, one needs to account 
for the free-space lost between the tangent plane and the 
sphere centred at the telescope and touching the plane 
(as shown in Figure]^. Both effects are negligibly small 
with the current experimental configuration. 

After calibrating all 32 tracks as described above, we 
have now a 2D plane with information about the beam on 
the grid formed by the tracks. We interpolate this plane 
using the python function scipy. interpolate. Rbf and 
epsilon=l to form a 2D map of the beam. The map is 
done for different frequency bins. 


4.2. ID beam pattern 

The ID beam pattern of the two tracks with maximum 
intensity is shown in Figure Each panel is constructed 
from an average of 20 frequency channels, with the mean 
frequency listed in the figure (1206 MHz, 1127 MHz and 
1012 M Hz, respect ively). To guide the eye, the Airy 


pattern (Airy 1838|) expected for an idealised 5m aper¬ 
ture with uniform illumination is overlaid in each panel 
with the black dashed curve. The thickness of the mea¬ 
surement curve corresponds to a 0.5 m error in the GPS 
positioning. 

First, it is worth noting that the measurements are at 
very high signal-to-noise even at the edge of the mea¬ 
surement (4th side-lobe). In principle, one can measure 


® These are typically regions where the interference from the 
drone actually raises the “off” signal to a non-negligible level. 


the beam pattern to further out with the current setup. 
We find that the Airy pattern gives a main beam size 
smaller than the drone measurement. This is expected 
as any de-focusing and aberration problems caused by 
imperfect geometry of the feed horn and the dish tend 
to enlarge the main beam. The nulls are in general not 
as deep as that predicted from the Airy pattern, but the 
positions of the first null and first side lobe agrees quite 
well. The positions of the higher-order side lobes and 
nulls in the measurement are shifted towards the main 
beam compared to the Airy pattern. The two measure¬ 
ments from the NS and EW tracks agree fairly well in 
the position and level of the peak/nulls, while the mea¬ 
surements show that the beam is not entirely symmetric. 
The asymmetry can be due to the intrinsic asymmetry of 
the telescope structure - the front-end unit is mounted 
slightly off-axis, and the dipole inside the feed horn is 
not centred. Finally, the beam size increases going to 
longe r wavelengths, as expected. We discuss further in 
§4.5| these wavelength-dependent characteristics. 

In Appendix]^ we invoke a simple antenna modelling 
tool GRASP to investigate the effect on the beam shape 
when the beam of the feed horn is included, using the 
geometry of the receiver horn described in §2.1| We also 
look at the impact on the beam shape from changes of 
other model parameters of the telescope. 

4.3. 2D beam pattern 

The reconstructed 2D beam pattern is shown in Fig- 
urej^for the same frequency ranges as Figure]^ Visually, 
one can see that the main beam is centred and roughly 
circular. The side-lobes show up as concentric ring struc¬ 
tures centred on the peak of the main beam, with dark 
rings indicating the nulls. There is noticeable asymme¬ 
try along the 45° direction. This is likely due to the po¬ 
larisation angle of the telescope (along the 45° NE-SW 
direction as shown in Figure H), breaking the otherwise 
isotropic beam pattern. TheM supporting struts could 
also introduce some of the structures in the beam pat¬ 
tern, as tehy are either parallel or perpendicular to this 
direction. Given the flight pattern design we used in this 
experiment (shown in Figure j^, artefacts from the grid- 
pattern are inevitable. That is, we only have data taken 
in certain stripes in the 2D plane, causing the reconstruc¬ 
tion to be limited in between the stripes, even if there 
were fine structures in the beam pattern. A more sophis¬ 
ticated flight pattern with adaptive grid size adjusted to 
the expected beam structure can potentially solve this 
issue in the future. 

These 2D beam maps can be generated for arbitrary 
frequency ranges and used as input to realistic simula¬ 
tions of the sky observed by the telescope. 

4.4. Other uncertainties 

We discuss and quantify here the possible sources of 
uncertainties in our results that we have not considered 
above. All of these effects are subdominant to that com¬ 
ing from the GPS uncertainty in the transverse direction. 

• Gimbal position/angle uncertainty: We mea¬ 
sured the stability of the angle of the gimbal to be 
within 1 deg, this corresponds to a <1 deg uncer¬ 
tainty in the polarisation of the emission from the 
transmitter, resulting in <0.1% uncertainty in the 
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Fig. 4. — Schematics of the flight pattern design (left) and the actual flight path during the experiment recorded by GPS (right). In the 
left panel, the black circle in the middle indicates the telescope, with the arrow showing the polarisation direction of the feed horn on the 
telescope. The blue (red) lines show the flight path in the NS (EW) orientation where signal from the noise transmitter is emitted. The 
dashed lines show the flight path where signal is not emitted. The triangle at the beginning of each dashed line shows the beginning of each 
flight, which include two tracks. The green circles and the arrows within indicates the drone and the orientation of the noise transmitter 
polarisation in the NS (lower green circle) and EW (upper green circle) direction. In the right panel, each stroke indicates an emission from 
the noise transmitter. The blue strokes are in the NS direction while the red strokes are the in EW direction.The black dots indicate the 
beginning of each track, corresponding to the triangles in the left panel. One can see that the distance between each stroke is not constant 
due to the change in the drone’s velocity. 
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Fig. 5.— Example of the data processing for the NS track closest to the beam centre. These curves correspond to a stacking of the 
frequencies in the range 1191.9 - 1220.2 MHz (20 channels with mean frequency 1206 MHz). Top row from left to right shows the raw 
data, the data adjusted for the units in both axes, the data corrected for the standing-wave pattern. The green arrows on the top left 
panel indicates the artefacts induced by the transmitter signal. Bottom row from left to right shows the data masking out the ’on’ and 
’off’ signal separately, interpolation over the data gaps, and the final corrected beam profile. 
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1206 MHz 






Fig. 6.— The centre beam profile measured in the NS and EW directions in linear (left) and log (right) scales for different frequencies. 
Each measurement curve is an average of 20 frequency channels, with the mean frequency being 1206 MHz (top), 1127 MHz (middle) and 
1012 MHz (bottom) respectively. Airy patterns corresponding to an idealised 5m aperture at the mean frequencies are shown by the black 
dashed curves. The width of the blue and red line indicate expected uncertainties from the imperfect GPS positioning. 




























1206 MHz 
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Fig. 7.— Reconstructed 2D beam intensity pattern from drone 
measurements for different frequencies. Each map is an average of 
20 frequency channels, with the mean frequency being 1206 MHz 
(top), 1127 MHz (middle) and 1012 MHz (bottom) respectively. 
The figure has the same orientation as Figure]^ Colour scales are 
logarithmic (dB). 

flux. The flux attenuation from the beam not di- 
rectly pointing at the dish corresponds to a <0.05% 
uncertainty in the flux level. 

• Radio Frequency Interference (RFI) re¬ 
moval: No severe high-level RFI contamination 
was observed during the flight, while low-level, long 
time-scale RFI is removed during the final step de¬ 
scribed in §4.1| 

• Polarisation uncertainty: One possibile error 
would arise if the flight path of the drone was not 


precisely 45 deg from the polarisation angle of the 
telescope. This would appear as a difference in the 
measured peak of the beam measured from the two 
main NS and EW tracks. We find this difference 
to be <4%. Note, however, that this number also 
contains the uncertainties in the background sub¬ 
traction, and other effects mentioned above. 

• Beam shape of transmitter horn: If the beam 
size of the transmitter horn is too small, it can im¬ 
print onto the beam measured from the telescope. 
We esti mate this effe ct using a simple horn antenna 
We find that for the tracks 


model (Kraus 


1950). 


furthest away from the dish centre, the flux atten¬ 
uation is about 8%, for the centre region where the 
main beam is probed (10 meters from the centre of 
the dish), the attenuation is below 1%. This means 
that we are underestimating the high-order slide 
lobes slightly, but the effect of the overall beam 
size measurement is small. 


4.5. Wavelength dependence 

In this section we calculate the wavelength-dependent 
beam characteristics from the measured beam profile. 
We take the two centre tracks used in Figure and keep 
all the frequency channels separate. Twenty lowest fre¬ 
quency channels were discarded due to severe RFI. The 
remaining frequency range is 997 - 1256 MHz. 

The beam FWHM (Equation is estimated by the 
FWHM of the best-fit Gaussian of the ID linear nor¬ 
malised profile. Hyi (Equation]^ and Dm (Equation]^ 
are calculated by integrating under the normalised beam 
map. In reality, we can only integrate Equation in¬ 
side our map (the inner 15 deg area of the beam). To 
account for the un-measured regions outside our map, 
we estimate the fraction of the beam inside our map 
from the beam model described in Appendix We 
find that ~ 85% of the beam is inside the 15 deg area. 
We thus multiply the measured integration of our maps 
1/0.85=1.18 to yield the used in Figurej^ As the first 
null in the single-frequency maps are often too noisy to 
define, for Dm we use the ID FWHM measurement to 
estimate the location of the first null in the case of an 
Airy pattern. From Figure!^ we can see that the position 
of the first null is quite wdl predicted by this approxi¬ 
mation. Finally, tja and tim can be calculated according 
to Equation and Equation 

The measured FWHM in the NS and EW directions 
as a function of frequency is shown in the first panel of 
Figure Also overlaid are the linear fits to both sets 
of data points. The fitted slopes are 1.28 and 1.25 for 
the NS and EW tracks respectively, which is about 20% 
larger than the idealised case (Equation . The plot 
also shows a prominent modulation alongrne linear re¬ 
lation. This modulation matches the time-independent 
standing-wave pattern observed in the data even without 
any signal present. These standing waves could be intro¬ 
duced b y multiple reflections off certain structures in the 
system (Briggs et al.||1997t [Popping & Braun||2008|. In 
our case, this is a combination of reflections from phys¬ 
ical components (e.g. dish-h orn) and reflections withi n 
the electronics. As shown in Popping & Braun (2008), 
the standing-wave imprints through the beam-size mea¬ 
surement, which is what we observe in the data. 
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Fig. 8. — Frequency-dependence of the beam characteristics. 
From top to bottom we plot the FWHM, beam efficiency and aper¬ 
ture efficiency as a function of frequency. The FWHM measure¬ 
ment is done on the ID profiles of the NS (blue) and EW (red) main 
track, while rfM and tja are calculated from the full 2D maps. 


The aperture efficiency rjM and the beam efficiency rjA 
as a function of frequency are shown in the second and 
third panel of Figure The median r]M over this fre¬ 
quency range is 68% and the median via is 67%. jja 
is fairly constant over frequency, while rjM increases by 
about 10% in our frequency range. Both plots also show 
imprints of the standing wave pattern, which appears to 
be out-of-phase with the FWHM measurement. This is 
expected, as a larger FWHM would lead to a smaller 
aperture/beam efficiency. 


4.6. Comparison with other measurements 

In this section we compare the drone measurements 
with a more traditional approach carried out at a time 
closely following the drone measurement, assuming that 
the beam shape is stable over that time period. It is 
also important that the instrument settings were kept 
the same as that used in the previous case. In §4.6.1| we 
repeat the measurement using the sun as a calibration 
source. We the n discu ss the pros and cons of the different 
approaches in §4.6.2| 


4.6.1. Beam measurement with the sun 

We performed the following sun scan on December 19, 
2014. The scanning strategy is designed to be similar to 
that of the drone measurement and allows us to recon¬ 
struct the beam pattern in both NS and EW directions. 

The data is taken with the telescope pointing South 
(azimuth 180°) and constantly changing elevation up and 
down from ^7° to ^32° elevation, the sun passes at el¬ 
evation around 20° while going east to west through az¬ 
imuth 180°. There are about 12 encounters of the sun 
and the beam where we see visible peak in the data. 
Zooming in each peak, one can see a full smooth profile 
from scanning the beam in the vertical direction, where 
the angular difference between the pixels is determined 
by the speed of the telescope slew. The peak of the pro¬ 
files can be identified as the point where the sun passed 
through the centre of that beam profile. Connecting all 
the peaks thus gives us the centre beam profile in the EW 
direction. Figure]^ shows the sun scan results of the two 
cross sections through the beam centre at a given fre¬ 
quency compared to the drone measurements. Results 
from other frequencies are similar. 

From Figurej^ we find that both the NS and the EW- 
beam from the sun are broadly consistent with the drone 
measurements at the l-2cr level. The SNR of the sun 
measurement is too low to resolve the side lobes except 
some hint of the first side lobe in both beams. 

We note that the main beam size from the sun mea¬ 
surement systematically smaller than that from the 
drone measurement. A few factors could contribute to 
this discrepancy: Eirst, the telescope had an elevation 
angle of 90° in the drone measurement and 19° in the sun 
measurement. This means that the mechanical structure 
could differ due to gravity and the level of ground pickup 
will be larger in the drone measurements. In addition, 
the drone signal was 10 dB larger than the sun signal, 
which suggested that any non-linear response from the 
instrument may also cause the two beams to produce 
different measurements. Quantifying exactly how much 
these contributes to the difference would however require 
more data taken over a longer period of time. 

4.6.2. Comparison of different measurement methods 

We have shown above that the two measurements of 
the beam are consistent, confirming that there are no un¬ 
known conceptual issues in using the drone to calibrate 
radio telescope beams. However, there are some funda¬ 
mental differences between the two different approaches 
of beam calibration. We summarise in Table the pros 
and cons of the two methods described above, together 
with other known techniques not covered in this paper. 
Overall, the drone measurement provides a more con¬ 
trollable way to calibrate the beam, and has potential 
to perform a broader range of calibrations (radiomet¬ 
ric calibration, beam as a function of elevation/azimuth, 
polarisation measurements etc.). On the other hand, the 
drone measurement becomes more challenging going to 
larger telescopes where the far-held is much further away 
from the antenna. Possible workarounds of this issue in¬ 
clude measuring the beam with the telescope pointing at 
a smaller elevation angle, or simply probing the near-held 
beam instead. 
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Fig. 9.— Linear (left) and log (right) beam profiles in the NS and EW cross sections through the beam centre as measured by the drone 
and the sun scan. The measurements correspond to an average of 20 frequency channels with the mean frequency 1206 MHz (same as the 
top panel In Figure]^. For the sun scan, the NS direction has a much lower sample rate due to the scanning strategy. 


TABLE 2 

Comparison between different approaches of radio beam calibration. 


Characteristic 

Drone 

Sun 

Weaker astronomical sources"' 

Available satellite^ 

Controllable flux, position and time 

yes 

no 

no 

no 

Radiometric calibration 

yes 

limited*^ 

yes 

yes 

Point source 

yes 

telescope-dependent*^ 

yes 

yes 

At infinity 

no 

yes 

yes 

yes 

SNR 

controllable 

lower 

lowest 

high 

Wavelength range 

broad 

broad 

depend on source 

limited 

Free 

no 

yes 

yes 

yes 


For example, moon, Cassiopeia A, Taurus A, Cygnus A, and Virgo A. 

Locations of the satellites in the sky need to be known beforehand. 

^ The flux and size of the sun varies with time, thus a good model of the sun is needed for accurately accounting its size and for 
radiometric calibration. 

Depending on the resolution of the telescope, the sun can be resolved in some cases. 


5. CONCLUSIONS 

In this paper, we describe a novel technique of cali¬ 
brating the beam of radio telescopes using drones. The 
advantage of this approach is that the calibration is con¬ 
trollable and flexible, which can be customised according 
to the focus of ones science goal. We demonstrate the 
approach by calibrating a 5m single dish at the Bleien 
Observatory in the 21 cm frequency range. We obtain 
high quality calibration data that allow us to understand 
the shape of the beam pattern in detail out to the 4th 
side-lobe. We characterise the wavelength dependency 
of the beam size, beam efficiency and aperture efficiency. 
The measurements are compared with more conventional 
measurements using the sun and the results are broadly 
consistent within measurement errors. We discussed the 
challenges in this experiment. Future improvements to 
the current experiment include: 

• Drone flight pattern design: The equal-spaced- 
grid flight pattern used in this work is not optimal, 
as there is more information in the centre of the 
beam. Using for example, an adaptive grid with 
varying spacing would be more effective in mapping 
the beam pattern. 

• Drone positioning system: One main error in 
our measurement comes from the inaccuracy of the 


GPS positioning system. This can be improved by 
using publicly available GPS augmentation systems 
or other techniques. 

• Characterisation of telescope: Our under¬ 
standing of the telescope geometry in this work 
is based on old drawings that could be outdated. 
Measuring more precisely the geometry of the tele¬ 
scope’s mechanical structure is important for bet¬ 
ter modelling. 

• Characterisation of horn feed: Measuring the 
beam of the horn alone in the lab would help dis¬ 
entangle the beam of the telescope from the total 
beam measured. This would also help us under¬ 
stand how the horn is illuminating the dish, and 
would be essential for further radiometric calibra¬ 
tions. 

• Modelling: In Appendix we use simple an¬ 
tenna modelling to help understand qualitative fea¬ 
tures in our measurements. However, in order to 
make quantitative comparisons between measure¬ 
ment and model, we require a more advanced soft¬ 
ware package. 

One of the science drivers for developing new, con¬ 
trollable beam calibration techniques for single-dish tele¬ 
scopes comes from the stringent requirement on the 
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knowledge of the telescope beam for cosmological HI in¬ 
tensity mapping. Small uncertainties in the beam would 
introduce undesirable systematics in the cosmological 
measurements. This work provides a practical solution to 
the challenge by building a controllable artificial calibra¬ 
tion source. This is achieved by combining commercial 
drone technology, well designed experiment setup and 
careful post processing of the data. 
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APPENDIX 


ANTENNA MODELLING 


In this appendix we use the software package GRASljjto produce a simple model of the expected beam pattern. GRASP 
is a standard engineering software used in designing reflector antennas. The main tuneable parameters we use in the 
software are the EM wave frequency, reflector geometry (surface type, diameter, f/D, offset) and horn feed geometry 
(taper angle, taper, polarisation). We are not able to add mechanical structures such as the supporting struts. 

The geometric information of the Bleien 5m dish can be obtained from the drawing of the telescope from the time of 
construction in the 1970’s, but this information does not include modifications on the telescope that have been done 
over the years. Uncertainty in the horn feed geometry is also present, as the current setup consists only of a low-cost 
cylindrical horn with a wire receiver mounted within. The horn was not designed to precisely match the beam of the 
telescope. We choose the default instrument parameters for our antenna model to be those listed in Table [ATT] These 
ar e set accordin g to the telescope drawing described above and an approximate model of the cylindrical horn described 


Silver ( 1984[). In our case, the dish is over-illuminated by the horn. 


We use GRASP mainly to understand the effect on the beam pattern when a more realistic bean of the feed horn is 
included. We also explore the impact of changing different parameters. Figure A.l for example, shows the effect of 
the b eam s hape when we perturb the tapering angle, the defocusing and the dish diameter from the fiducial setting in 


Table A.l We find that change in these parameters have an impact on the positions of the peaks/n ulls as well as the 


relative height of the different peaks. These qualitative changes are consistent with that described in Baars (20071. A 
more sophisticated model is needed to address more subtle changes in the shape of the beam. 


® http://www.ticra.com/products/software/grasp 
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TABLE A.l 

Default parameters for 

GRASP ANTENNA MODELLING. 


Parameter 

Value 

Reflector type 

parabolic 

Reflector f/D 

0.507 

Reflector diameter 

5m 

Frequency 

1114 MHz 

Feed taper angle 

150° 

Feed taper 

-10 dB 

Defocus 

0 
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Fig. A.l.— Beam patterns from GRASP when various instr umen t parameters are varied. The black dashed curve in all panels are the 
same, and show the fiducial beam pattern according to Table [AT| The three rows show the effect of varying defocus, tapering angle and 
the dish size, respectively. The four numbers listed in the legend are [taper angle, taper, radius (m), defocus (cm)] and describe the change 
in parameters. 




























































